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IN THE SPECIFICATION: 

Please amend the specification as follows: 

In the paragraph beginning at page 14, line 24: 

Other objects and advantages of the present invention will become apparent from the 
following description taken in conjunction with the accompanying drawings, which constitute a 
part of this specification disclosure and wherein are set forth exemplary embodiments of the 
present invention to illustrate various objects and features thereof 

In the paragraph beginning at page 10. line 13: 

(f) Since the cubic splines can have wide swings at the ends, the envelope- 
mean method of the EMD method is particularly unsuitable for real-time applications or for 
applications utilizing a narrow window. The end effects also propagate to the interior and 
significantly corrupt the data, as the construction of MFs progresses, as can be seen in Figures 1- 
2. Figure 1, left panel, illustrates the application of the EMD to a test signal (top-most signal) to 
produce IMF components (displayed below test signal). This panel illustrates end effects, spline- 
related instabilities (most noticeable in bottom components), and inability to extract the readily 
apparent signal baseline. The intrinsic timescale decomposition, sometimes referred to herein as 
ITD 4 (shown in the right panel) separates the same test signal (top-most signal) into stable 
components (displayed below test signal), demonstrating the fact that that ITD has no end effect 
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propagation beyond the first two extrema at each level and allows correct identification of the 
trend (dashed line). Figure 2, top panel, shows a brain wave input signal (electrocorticogram; 
'ECoG') containing an epileptic seizure used to illustrate decomposition differences between 
EMD and ITD. EMD-generated (lower left panel) and ITD-generated (lower right panel) 
decompositions of the cumulative sum of the raw signal show that ITD, unlike EMD, does not 
generate extraneous components and correctly reveals large timescale variations of the signal 
(DC trend). 

In the paragraph beginning at page 32. line 7: 

Figure 5 demonstrates how the ITD may be applied to determine a TFE distribution of a 
sample signal and compares its performance to prior art methods of Fourier analysis and wavelet 
analysis for this task. Figure 5A contains a sample signal. Figures 5B-C contain the TFE 
distribution of the signal obtained via windowed Fourier analysis. Note the rectangular grid of 
TFE information naturally derived from this transform and the trade-offs between good temporal 
but poor frequency localization (as shown in Figure 5B), and poor temporal but good frequency 
localization (as shown in Figure 5C). Darker regions correspond to those with higher power 
throughout this Figure. It is impossible to simultaneously improve both temporal and frequency 
localization of energy. Quantification is performed in this Fourier analysis by sequential time 
windows with a predetermined segmentation of the signal that is not determined in any way by 
the signal changes (e.g., local extrema). Figure 5D illustrates the TFE distribution for the same 
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input signal that is obtained using the fast wavelet transfor m, sometimes referred to herein as the 
FWT . Note the dyadic grid of TFE information naturally derived from this transform. The FWT 
has the ability to localize higher frequency information on shorter time scales and lower 
frequency information on longer timescales, but still suffers from predetermined temporal 
segmentation and inaccuracies in both temporal and frequency localization. These are due in part 
to the bleeding of signal energy across different levels of resolution in the wavelet transform and 
are attributed in large part to the temporal and frequency bins that are predetermined by the 
choice of basis and not by the signal under study. Figure 5E shows the ITD-based TFE 
distribution for the same signal, illustrating the significant advance in time-frequency localization 
provided by the algorithm of the subject invention. Each displayed line segment's start and end 
points correspond exactly to the start and end of an actual wave in an ITD proper rotation 
component and the line segment shading is determined by the amplitude of the corresponding 
wave (darker indicates greater amplitude). Since the ITD automatically decomposes the input 
signal into a set of proper rotations and a monotonic trend, the instantaneous phase angles for 
each proper rotation component can be obtained via the Hilbert transform or the invention's 
alternative approaches described herein to obtain the corresponding instantaneous frequencies. 
Figure 5F (page 2 of 3) shows the original signal and the first four HF components obtained by 
the ITD. Figure 5G shows the instantaneous frequency curves (in light gray when instantaneous 
power is insignificant, i.e., below 2xl0" 5 ) corresponding to each of the HF proper rotation 
components (as shown smoothed and solid black) when instantaneous power is non-negligible, 
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i.e., above 0.005, to illustrate this process. The curves in Figure 5G were obtained using the 
arcsine-based method for instantaneous frequency computation and were smoothed using a 0.1 
second median filter to enhance visualization. Figure 5H (page 3 of 3) shows the first three 
proper rotation components obtained by application of the ITD to the signal of Figure 5 A. Figure 
51 and 5J show the half-wave amplitudes and (smoothed) instantaneous frequencies for each 
component obtained via the arcsine-based method for instantaneous phase and instantaneous 
frequency determination (linestyles are the same as that of the corresponding proper rotation 
components). As expected, the instantaneous frequencies of each component decrease at each 
point in time as the level increases, but the actual frequency values are determined solely by the 
signal and not predetermined in any manner. 

In the description of Fig. 1 on page 15: 

Figure. 1 Figures la- lb is a comparison of the prior art method of EMD and with the ITD 
as applied to a test signal. 



In the description of Fig. 2 on page 1 5: 

Figured Figures 2b-2c is a comparison of the prior art method of EMD and with the ITD 
as applied to a brain wave signal as shown in Fig . 2a 
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In the description of Fig. 4 on page 1 5: 

Figure 4 is an illustiatiun Figures 4a-4 d are illustrations of the steps of the online ITD 
method. 

In the description of Fig. 5 on page 1 5: 

Figure 5 is a eompaiisun Figures 5b-5d are comp arisons between the ITD and prior art 
methods of Fourier analysis and wavelet analysis and the ITD TFjgs. 5e-5i^ in determining time- 
frequency-energy (TFE) distributions from a sample signal as shown in Fig. 5a . 

In the description of Fig. 6 on page 1 5: 

Figure 0 is an illuslialiun Figures 6a-6 g are illustrations of the ITD-based method. 

In the description of Fig. 7 on page 1 5: 

Figure 7 ii an illuslialiun es Figures 7a-7h are illustrations of the ability of the ITD-based 
filtering method, used to differentiate between two types of waves that have significantly 
overlapping spectral characteristics. 

In the description of Fig. 8 on page 1 5: 

Figure 0 is an illusliali u n Figures 8a-8c are illustrations of ITD applied to 
multidimensional data. 
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After the description o f Figure 8 on page 1 5. insert thefollowing : 

Figure 9A is a flow chart representation of the ITD method of the present invention 

Figure 9B is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9B of Figure 9 A. 

Figure 9C is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9C of Figure 9 A. 

Figure 9D is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9D of Figure 9A. 

Figure 9E is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9E of Figure 9 A. 

Figure 9F is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9F of Figure 9A. - - 

In the paragraph beginni ng at page 27. line 1 V 

Figure 4 illusudles Figures 4a-4d illustrate the procedure followed by the online ITD 
system. In Figure 4A 4a, the input signal to be decomposed is shown in solid lines with extrema 
indicated by circles. The baseline component, computed in time up to the local minimum 
extrema at time t(j-l) is shown in dashed lines. The system is buffering each new input signal 
value on the solid curve until such time as a new extrema is detected at time t(j+ 1). The 
detection of the new extrema triggers computation (via the ITD step algorithm) of the monotonic 
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baseline segment for values of t between t(j-l) and t{j). First, a new baseline node is computed 
(see solid point at t(j) in Figure 4B 4b), and the input signal itself over that interval is linearly 
transformed to form the monotonic baseline function over that interval (see dashed red curve 
between t(j-l) and t(j) in Figure 4€ 4c). The low frequency baseline component segment is then 
immediately subtracted from the input signal to generate the corresponding high frequency 
component on the \tQ-l), t(j)} interval. The procedure continues, with each new extrema 
triggering a similar portion of the decomposition on the time interval between adjacent extrema. 
See Figure 4© 4d for an illustration of resulting baseline and high frequency components for this 
example, along with the original signal being decomposed. Note that: (i) the high frequency 
component has all positive maxima and all negative minima (i.e., it is a "proper rotation"), and 
(ii) the method is iterated, applying this decomposition procedure to each resulting baseline 
component in a manner that is similarly triggered each time the component to be decomposed has 
a new extrema. 

In the para graph beginning at page 32. line 9: 

Figure 5 dunuiisUdtes Figures 5a-5i demonstrate how the ITD may be applied to 
determine a TFE distribution of a sample signal and compares its performance to prior art 
methods of Fourier analysis and wavelet analysis for this task. Figure 5A- 5a contains a sample 
signal. Figures 5B=6 5b-5c contain the TFE distribution of the signal obtained via windowed 
Fourier analysis. Note the rectangular grid of TFE information naturally derived from this 
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transform and the trade-offs between good temporal but poor frequency localization (as shown in 
Figure 5B 5b ), and poor temporal but good frequency localization (as shown in Figure 5€ 5c ). 
Darker regions correspond to those with higher power throughout this Figure. It is impossible to 
simultaneously improve both temporal and frequency localization of energy. Quantification is 
performed in this Fourier analysis by sequential time windows with a predetermined 
segmentation of the signal that is not determined in any way by the signal changes (e.g., local 
extrema). Figure 5© 5d illustrates the TFE distribution for the same input signal that is obtained 
using the fast wavelet transform. Note the dyadic grid of TFE information naturally derived from 
this transform. The FWT has the ability to localize higher frequency information on shorter time 
scales and lower frequency information on longer timescales, but still suffers from predetermined 
temporal segmentation and inaccuracies in both temporal and frequency localization. These are 
due in part to the bleeding of signal energy across different levels of resolution in the wavelet 
transform and are attributed in large part to the temporal and frequency bins that are 
predetermined by the choice of basis and not by the signal under study. Figure SE 5e shows the 
ITD-based TFE distribution for the same signal, illustrating the significant advance in time- 
frequency localization provided by the algorithm of the subject invention. Each displayed line 
segment's start and end points correspond exactly to the start and end of an actual wave in an 
ITD proper rotation component and the line segment shading is determined by the amplitude of 
the corresponding wave (darker indicates greater amplitude). Since the ITD automatically 
decomposes the input signal into a set of proper rotations and a monotonic trend, the 
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instantaneous phase angles for each proper rotation component can be obtained via the Hilbert 
transform or the invention's alternative approaches described herein to obtain the corresponding 
instantaneous frequencies. Figure 5F 5f (page 2 of 3) shows the original signal and the first four 
HF components obtained by the ITD. Figure 5G 5g shows the instantaneous frequency curves (in 
light gray when instantaneous power is insignificant, i.e., below 2x10 s ) corresponding to each of 
the HF proper rotation components (as shown smoothed and solid black) when instantaneous 
power is non-negligible, i.e., above 0.005, to illustrate this process. The curves in Figure 56 5g 
were obtained using the arcsine-based method for instantaneous frequency computation and were 
smoothed using a 0.1 second median filter to enhance visualization. Figure 5H 5h (page 3 of 3) 
shows the first three proper rotation components obtained by application of the ITD to the signal 
of Figure 5A 5a . Figure 51 Figures Si and 5J 5j show the half-wave amplitudes and (smoothed) 
instantaneous frequencies for each component obtained via the arcsine-based method for 
instantaneous phase and instantaneous frequency determination (linestyles are the same as that of 
the corresponding proper rotation components). As expected, the instantaneous frequencies of 
each component decrease at each point in time as the level increases, but the actual frequency 
values are determined solely by the signal and not predetermined in any manner. 



In the paragraph begin ning at page 34. line 6: 

The ability of the ITD to decompose online any input signal into a set of one or more 
component signals which have the proper rotation property makes it an ideal first step of a two- 
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step process resulting in a powerful and unique new method for nonlinear signal filtering. In 
particular, as each monotonic segment is computed within one or more decomposition levels, 
features of the segment can be analyzed and quantified. For example, the amplitude and duration 
of the segment may be simply computed from the extrema, i.e., the starting and ending points of 
the segment. One skilled in the art will appreciate that more detailed features such as the time- 
weighted feature density of the segment, see for example Int'l. Publication No. WO 01/75660 
Al, of Nikitin et al, and associated properties such as its average or median value, its variance, 
skewness, distance from a template or reference density, etc., may also be quantified for each 
segment. Individual segments may then be classified according to the values of theses quantified 
features. A filtered signal can then be constructed by adding together (i.e., superimposing the 
segments while preserving the time intervals on which each occurred) only those segments that 
satisfy certain constraints applied to their quantified features. For example, one may construct a 
signal using only those ITD-created monotonic segments with durations that, using the data 
sampling rate, correspond to wavespeeds between f^ Hz and f^ Hz, and with absolute 
amplitudes of at least A^, e.g., in excess of the 75 th percentile of amplitudes for all segments 
with the required duration. In this example, the output can be essentially interpreted as extracting 
the waves in the input signal that are in the f^ - f^ frequency band and which have larger 
amplitudes than A min . Such uses would bear a strong resemblance to the commonly used Fourier- 
based band-pass filtering, but would avoid the phase-shifting and waveform distortion 
drawbacks, while adding the amplitude differentiation capability. In Figur e 6 Figures 6a-6g . this 
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concept is illustrated when applied to the same input signal used in Figure SA 5a (also shown 
again in Figure 6A 6a ). After the TFE distribution is computed, the absolute amplitude and 
duration/wavespeed of each individual wave are obtained. Figur e , 6D-C shows Figures 6b-6c 
show the one-dimensional (marginal) densities of these two features along with vertical 
discriminator in solid lines that were selected for purposes of illustrating the capabilities of this 
method. Figure 6B 6d shows the (absolute amplitude, wavespeed) ordered pairs for each 
segment. Discriminators of pairs with wavespeed less than 20 Hz, between 20 Hz and 50 Hz, 
and greater than 50 Hz are applied, along with a discriminator of amplitude exceeding 0.01 mV. 
Feature vector pairs (i.e., points in Figure 6© 6d ) falling into certain specific regions of the two- 
dimensional feature range are classified according to the region in which they belong. Figure 6E 
6e displays three filtered components constructed using only those segments with points in 
specific regions (corresponding to three of the regions in Figure 6© 6d ). Component one, C„ is 
the output of a filter that contains the signal reconstructed using only those waves that have 
absolute amplitude exceeding O.OlmV and wavespeed between 20Hz and 50 Hz. Component 
two, C 2 , is the output of a filter that contains the signal reconstructed using only those waves that 
have absolute amplitude less than O.OlmV and wavespeed between 20Hz and 50 Hz. Component 
three, C 3 , is the output of a filter that contains the signal reconstructed using only those waves 
that have absolute amplitude exceeding O.OlmV and wavespeed below 20Hz. The output 
obtained by combining some of these components (i.e., C,+C 2 +C 3 and C,+C 2 ) are shown in 
Figure 6F 6f , along with the output of a traditional, linear, [20, 50] Hz. band-pass digital filter 
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provided for comparison. The extrema-preserving properties of this new type of filter are evident 
in this Figure 6F 6f. Figure 6G 6g again compares an output available by using the ITD filter in 
comparison to the same linear band-pass digital filter, but does so on a slightly larger timescale 
(10 seconds). The residual signals obtained by subtracting the filtered outputs from the original 
signal are also shown, to illustrate the ability of the ITD based signal to more completely 
decompose the raw input signal into what may be interpreted as abnormal and normal 
components present in the raw signal. 



In the par agraph beginning at page 36. line 24: 

Another example that illustrates the power of the invention for signal filtering and 
analysis is provided in Figttrer? Figures 7a-7h. In thuflguu, these Figures, we demonstrate how 
the ITD-based filtering method may be applied to differentiate between two types of waves that 
have significantly overlapping spectral characteristics and to decompose a raw signal that is a 
combination of these two types of waves into components that preserve temporal location of 
extrema and critical points in the raw signal. Overlap between the power spectral densities of an 
underlying signal and superimposed noise is one of the more difficult problems in signal analysis 
and detection, with numerous applications. Figure 9A 7a illustrates the mix of two alternating 
signals, one a 9 Hz cosine wave with an increasing linear trend and the other a 9 Hz sawtooth 
wave. The signals were added to create a raw test signal for use in this example. Figure 9B 7b 
shows the power spectral density estimates for each of the two signals, clearly illustrating the 
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significant overlap in the band around 9 Hz. Figure 9B 7c shows the proper rotation components 
and monotonic trend obtained via application of ITD to the original combined signal. Figure 7D- 
E Figures 7d-7e illustrate the result of applying single wave analysis to examine various features 
of individual waves obtained by the ITD decomposition. Figure m 7d shows the ratio of 
individual wave peak to mean signal value, while Figure 9E 7e shows a plot of wave skewness 
plotted versus wave kurtosis. A sample discriminator in dashed linestyle is also shown to 
illustrate the result of feature classification as applied to these quantified waveform features. 
Figure 7F 7f shows the output of the ITD-based filter that retains only those waves with peak-to- 
mean value ratio exceeding 1.65. The corresponding residual signal, i.e., the original signal 
minus the filter output signal, is also shown and illustrates that peak-to-mean value ratio 
discrimination of ITD-produced proper rotation waves allow reconstruction of a filter output that 
nearly perfectly decomposes the original signal into its constituent components while retaining 
precise temporal information regarding, e.g., extrema location. Figure 9B 7g shows the output 
of the ITD-based filter that retains only those waves with skewness-kurtosis pairs above the 
dashed discriminator line, along with the corresponding residual signal. This example illustrates 
how this new filter method may be easily configured to decompose an input signal into 
components with more complex criteria, such as "cutting off the tops of all riding sinusoidal 
waveforms while retaining every other signal component" as is done here. Figure m 7h shows 
the same signals as in Figure 7G but zooming in to the time interval between t=9 and t=13 to see 
details of the filter output. 



Application No. 10/684,189 
Mark G. Frei et al, Applicants 
Page 18 

In the par agraph beginning at page 39. line 6: 

One skilled in the art will also recognize that, just as Fourier-based filtering and other 
decomposition methods can be successfully applied to higher dimensional analysis, so, too, can 
the ITD system. This extends the above-mentioned application of the ITD to decompose 
functions of several variables. For example, in image processing applications, one may be 
interested in separating an image into a high frequency component, which contains edge 
transitions between various objects in the image for example, and a low frequency component, 
which define background colors in the image for example. Two-dimensional wavelet and/or 
Fourier analysis are popular for image processing applications. However, the drawbacks 
mentioned above for these methods still exist in higher dimensional analysis. The benefits of the 
ITD system, such as its ability to preserve extrema and precisely localize time-frequency-energy 
information are preserved in higher dimensional analysis and thereby offer improvements over 
the most popular prior art methods. Figure 0 ihWiatu, Figures 8a-8c illustrate the application of 
ITD to decompose a two-dimensional surface (shown in Figure 8A 8a) into two component 
surfaces, namely a high frequency surface (shown in Figure 8B 8b) and a low frequency surface 
(shown in Figure 86 8c). The decomposition in this example is obtained by decomposing each 
of the cross-sectional signals of the surface along the grid lines obtained by holding the first 
independent variable constant, repeating the process holding the second variable constant, and 
averaging the two results. 
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In the paragraph beg i nning at page 34. line 6: 

The ability of the ITD to decompose online any input signal into a set of one or more 
component signals which have the proper rotation property makes it an ideal first step of a two- 
step process resulting in a powerful and unique new method for nonlinear signal filtering. In 
particular, as each monotonic segment is computed within one or more decomposition levels, 
features of the segment can be analyzed and quantified. For example, the amplitude and duration 
of the segment may be simply computed from the extrema, i.e., the starting and ending points of 
the segment. One skilled in the art will appreciate that more detailed features such as the time- 
weighted feature density of the segment , see for example Int'l. Publication Nu. WO 01/75660 
Al, of Nikilin ct al, and associated properties such as its average or median value, its variance, 
skewness, distance from a template or reference density, etc., may also be quantified for each 
segment. Individual segments may then be classified according to the values of theses quantified 
features. A filtered signal can then be constructed by adding together (i.e., superimposing the 
segments while preserving the time intervals on which each occurred) only those segments that 
satisfy certain constraints applied to their quantified features. For example, one may construct a 
signal using only those ITD-created monotonic segments with durations that, using the data 
sampling rate, correspond to wavespeeds between Hz and f^ Hz, and with absolute 
amplitudes of at least A,^, e.g., in excess of the 75* percentile of amplitudes for all segments 
with the required duration. In this example, the output can be essentially interpreted as extracting 
the waves in the input signal that are in the f min - f max frequency band and which have larger 
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amplitudes than A^^. Such uses would bear a strong resemblance to the commonly used Fourier- 
based band-pass filtering, but would avoid the phase-shifting and waveform distortion 
drawbacks, while adding the amplitude differentiation capability. In Figure 6, this concept is 
illustrated when applied to the same input signal used in Figure 5 A (also shown again in Figure 
6A). After the TFE distribution is computed, the absolute amplitude and duration/wavespeed of 
each individual wave are obtained. Figure 6B-C shows the one-dimensional (marginal) densities 
of these two features along with vertical discriminator in solid lines that were selected for 
purposes of illustrating the capabilities of this method. Figure 6D shows the (absolute amplitude, 
wavespeed) ordered pairs for each segment. Discriminators of pairs with wavespeed less than 20 
Hz, between 20 Hz and 50 Hz, and greater than 50 Hz are applied, along with a discriminator of 
amplitude exceeding 0.01 mV. Feature vector pairs (i.e., points in Figure 6D) falling into certain 
specific regions of the two-dimensional feature range are classified according to the region in 
which they belong. Figure 6E displays three filtered components constructed using only those 
segments with points in specific regions (corresponding to three of the regions in Figure 6D). 
Component one, C„ is the output of a filter that contains the signal reconstructed using only 
those waves that have absolute amplitude exceeding O.OlmV and wavespeed between 20Hz and 
50 Hz. Component two, C 2 , is the output of a filter that contains the signal reconstructed using 
only those waves that have absolute amplitude less than O.OlmV and wavespeed between 20Hz 
and 50 Hz. Component three, C 3 , is the output of a filter that contains the signal reconstructed 
using only those waves that have absolute amplitude exceeding O.OlmV and wavespeed below 
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20Hz. The output obtained by combining some of these components (i.e., C,+C 2 +C 3 and C,+C 2 ) 
are shown in Figure 6F, along with the output of a traditional, linear, [20, 50] Hz. band-pass 
digital filter provided for comparison. The extrema-preserving properties of this new type of 
filter are evident in this Figure 6F. Figure 6G again compares an output available by using the 
ITD filter in comparison to the same linear band-pass digital filter, but does so on a slightly 
larger timescale (10 seconds). The residual signals obtained by subtracting the filtered outputs 
from the original signal are also shown, to illustrate the ability of the ITD based signal to more 
completely decompose the raw input signal into what may be interpreted as abnormal and normal 
components present in the raw signal. 



